Determining Cosmological Parameters with Latest Observational Data 
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In this paper, we combine the latest observational data, including the WMAP five-year data 
(WMAP5), BOOMERanG, CBI, VSA, ACBAR, as well as the Baryon Acoustic Oscillations (BAO) 
and Type la Supernoave (SN) "Union" compilation (307 sample), and use the Markov Chain Monte 
Carlo method to determine the cosmological parameters, such as the equation-of-state (EoS) of 
dark energy, the curvature of universe, the total neutrino mass and the parameters associated 
with the power spectrum of primordial fluctuations. Our results show that the ACDM model 
remains a good fit to the current data. In a flat universe, we obtain the tight limit on the constant 
EoS of dark energy as, w — —0.977 ± 0.056 (1 a). For the dynamical dark energy models with 
time evolving EoS parameterized as Wdc{a) = wo + ■wi{l — a) , we find that the best-fit values are 
Wo = —1.08 and wi — 0.368, implying the preference of Quintom model whose EoS gets across the 
cosmological constant boundary during evolution. For the curvature of universe Qk, our results give 
—0.012 < Qk < 0.009 (95% C.L.) when fixing lUdc = —1. When considering the dynamics of dark 
energy, the flat universe is still a good fit to the current data, —0.015 < flk < 0.018 (95% C.L.). 
Regarding the neutrino mass limit, we obtain the upper limits, ^ m^ < 0.533 eV (95% C.L.) within 
the framework of the flat ACDM model. When adding the SDSS Lyman-a forest power spectrum 
data, the constraint on ^ m^, can be significantly improved, ^ m^ < 0.161 eV (95% C.L.). However, 
these limits can be weakened by a factor of 2 in the framework of dynamical dark energy models, 
due to the obvious degeneracy between neutrino mass and the EoS of dark energy model. Assuming 
that the primordial fluctuations are adiabatic with a power law spectrum, within the ACDM model, 
we find that the upper limit on the ratio of the tensor to scalar is r < 0.200 (95% C.L.) and the 
inflationary models with the slope n^ > 1 are excluded at more than 2 a confidence level. However, 
in the framework of dynamical dark energy models, the allowed region in the parameter space of 
{ns,r) is enlarged significantly. Finally, we find no strong evidence for the large running of the 
spectral index, as = -0.019 ± 0.017 (1 a) for the ACDM model and a^ = -0.023 ± 0.019 (1 a) for 
the dynamical dark energy model, respectively. 



I. INTRODUCTION 

With the accumulation of observational data from cos- 
mic microwave background measurements (CMB), large 
scale structure surveys (LSS) and supernovae observa- 
tions and the improvements of the data quality, the cos- 
mological observations play a crucial role in our under- 
standing of the universe, and also in constraining the 
cosmological parameters, such as the EoS of dark en- 
ergy models, the curvature of universe, the total neu- 
trino masses and those associated with the running of 
the spectral index and gravitational waves. In our pre- 
vious work [l|, we have used the Markov Chain Monte 
Carlo (MCMC) method to constrain cosmological mod- 
els from the astronomical observational data, including 
the WMAP three-year data (WMAP3) [3,[1, small-scale 
CMB data, LSS data [1 i] and SNIa ESSENCE sample 
[a|. We found that the cosmological constant is consis- 
tent with the data, however the dynamical dark energy 
models are still allowed and interestingly the model with 
its EoS getting across w = —1 the Quintom model 7] is 
the best fit model. And we found no strong significant 
evidence for the non-flat universe and massive neutrino. 
Within the ACDM model, the scale-invariant spectrum 
and the spectra with n^ > 1 are disfavored by more than 
2 a confidence level. Due to the degeneracy between the 



EoS of dark energy and tensor fluctuation, those infla- 
tionary models excluded within the ACDM model, will 
revive in the framework of dynamical dark energy mod- 
els. Furthermore, we did not find any significant evidence 
for the tensor fluctuations and the large running of the 
spectral index. 

Given the precision of current observations, these re- 
sults are not conclusive for the time being. Recently, 
the WMAP experiment has published its five-year data 
of temperature and polarization power spectra [^, |9|, [lOj . 
The Arcminute Cosmology Bolometer Array (ACBAR) 
experiment has also published its new CMB tempera- 
ture power spectrum [ll|. These new CMB data can 
strengthen the constraints on the cosmological parame- 
ters, especially for the inflationary models [1, l^jll4|- Fur- 
thermore, the Supernova Cosmology Project has made an 
unifled analysis of the world's supernovae datasets and 
presented a new compilation "Union" (307 sample) [1J| 
which includes the recent samples of SNIa from SNLS 
and ESSENCE Survey, as well as some older datasets, 
etc. In the literature [ij, [ij], this Union compilation 
combining with the shift parameter derived from CMB 
and the BAO information has been used to constrain cos- 
mological models. However, in these studies the CMB in- 
formation considered is just the shift parameter instead 
of the full CMB data, which will lose some information 



of CMB and lead to a biased result [l5|. Thus, in this 
paper, we revisit the issue on the determination of these 
cosmological parameters and update our previous results 
with the latest observational data. 

Our paper is organized as follows: In Section II we de- 
scribe the method and the latest observational datasets 
we used; Section III contains our main global fitting re- 
sults on the cosmological parameters and the last section 
is the summary. 



II. METHOD AND DATA 

In our study, we perform a global analysis using the 
publicly available MCMC package CosmoMC^ [ll|. We 
assume the purely adiabatic initial conditions. Our most 
general parameter space is: 

P = {u!b,uJc,^k,Qs,T,wo,wi, fu,ns,As,as,r) , (1) 

where uji, = fibft,^ and oJc = ^ch^, in which f2{, and flc 
are the physical baryon and cold dark matter densities 
relative to the critical density, ilfc is the spatial curva- 
ture and satisfies Qk + ^m + ^de = 1, Os is the ratio 
(multiplied by 100) of the sound horizon to the angular 
diameter distance at decoupling, r is the optical depth 
to re-ionization, f^ is the dark matter neutrino fraction 
at present, namely. 
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where A^ is defined as the amplitude of initial power 
spectrum, rig measures the spectral index, as is the run- 
ning of the scalar spectral index and r is the tensor to 
scalar ratio of the primordial spectrum. For the pivot 
scale we set kso = 0.05Mpc~-'^. Moreover, wq and wi are 
the parameters of dark energy EoS, which is given by 



Wdcia) ^ Wo + wi{l - a) 



(4) 



where a — 1/(1 + z) is the scale factor and wi — —dw/da 
characterizes the "running" of EoS (RunW henceforth). 
The ACDM model has wq = -1 and wi = 0. For the 
dark energy model with a constant EoS, wi = (WCDM 
henceforth). When using the global fitting strategy to 



constrain the cosmological parameters, it is crucial to 
include dark energy perturbations 0, U^ |20|. In this 
paper we use the method provided in Refs.|20l. l2ll | to 
treat the dark energy perturbations consistently in the 
whole parameter space in the numerical calculations. 

In the computation of CMB we have included the 
WMAP5 temperature and polarization power spectra 
with the routine for computing the likelihood supplied 
by the WMAP team^. We also include some small- 
scale CMB measurements, such as BOOMERanG |23], 
CBI [^, VSA m and the newly released ACBAR data 
[ll| . Besides the CMB information, we also combine 
the distance measurements from BAO and SNIa. For 
the BAO information, we use the gaussian priors on the 
distance ratios, rJDyiz) = 0.1980 ± 0.0058 at z = 0.2 
and rs/Dy{z) = 0.1094 ± 0.0033 at z = 0.35, with a 
correlation coefficient of 0.39, measured from the BAO 
in the distribution of galaxies [23|. In the calculation 
of the likelihood from SNIa we have marginalized over 
the nuisance parameter |26| . The supernova data we 
use are the recently released "Union" compilation (307 
sample) [13|. In order to improve the constraint on 
the total neutrino mass, we include the Lyman-a forest 
30wer spectrum from Sloan Digital Sky Survey (SDSS) 
27|, however also keep its unclear systematics in mind 
|28|. Furthermore, we make use of the Hubble Space 
Telescope (HST) measurement of the Hubble parame- 
ter Hq = lOOh km s^^ Mpc~^ by a Gaussian likelihood 
function centered around h — 0.72 and with a standard 
deviation CT = 0.08 [H. 



III. NUMERICAL RESULTS 

In this section we present our global fitting results of 
the cosmological parameters determined from the latest 
observational data and focus on the dark energy param- 
eters, curvature of universe, neutrino mass and the infia- 
tionary parameters, respectively. 



A. Equation of State of Dark Energy 

In Table I we list the constraints on the dark energy 
parameters as well as the Hubble constant in different 
dark energy models. 

Assuming the fiat universe, firstly we explore the con- 
straints on the constant EoS of dark energy, w (w = 
Wq, Wi =0), from the latest observational data. In Fig[T] 
we show the constraints on w and the present dark mat- 
ter density, Qm- This result shows that the combina- 
tion of Union compilation (with systematic uncertain- 
ties not included) and other observational data yield a 
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TABLE I. Constraints on the dark energy EoS and some background parameters from the latest observations. Here we have shown the 
mean and the best fit values, which are obtained from the cases with and without the systematic uncertainties of Union compilation, 

respectively. 



Parameter 


WQ 


Wl 


Ode 


Ho 




with sys 


w/o sys 


with sys 


w/o sys 


with sys 


w/o sys 


with sys 


w/o sys 


ACDM 


BestFit 


-1 


-1 








0.735 


0.741 


71.0 


71.6 


Mean 


-1 


-1 








0.738 ±0.015 


0.738 ±0.014 


71.4 ±1.4 


71.4 ±1.3 


WCDM 


BestFit 


-0.978 


-0.955 








0.738 


0.735 


71.4 


70.6 


Mean 


-0.965 ± 0.080 


-0.977 ±0.056 








0.736 ±0.016 


0.737 ±0.014 


70.8 ±1.9 


71.1 ±1.4 


RunW 


BestFit 


-1.09 


-1.08 


0.533 


0.368 


0.735 


0.738 


70.4 


71.1 


Mean 


-0.946 ±0.194 


-0.993 ±0.128 


-0.133 ±0.749 


0.030 ±0.582 


0.734 ±0.017 


0.737 ±0.014 


70.7 ±1.9 


70.9 ±1.5 


RunW 


BestFit 


- 


-1.11 


- 


0.475 


- 


0.739 


— 


72.4 


Mean 


- 


-0.976 ±0.148 


- 


-0.071 ±0.848 


- 


0.736 ±0.014 


- 


70.9 ±1.9 


RunW 
w/o Pert. 


BestFit 


- 


-1.04 


- 


0.290 


- 


0.742 


- 


71.3 


Mean 


- 


-1.00 ±0.114 


- 


0.103 ±0.413 


- 


0.736 ±0.012 


- 


70.8 ±1.5 
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FIG. 1: Constraints on the constant EoS of dark energy, w, 
and the present matter density, flm, from the latest observa- 
tions, assuming a flat universe. The blue dash-dot lines and 
red solid lines are obtained with and without the systematic 
uncertainties of Union compilation, respectively. 



strong constraint on the constant EoS of dark energy, 
w = -0.977 ± 0.056 (la). Our result is similar to the 
limit from WMAP5 [8], w = -0.972+°°^J (1 a) and 
physically it indicates that w = — 1 is consistent with 
the data. Furthermore, some of the quintessence mod- 
els get strongly constrained, for example, the tracker 
quintessence model which predicts w ~ —0.7 [30| will be 
excluded by more than 5 a. However we notice that the 
systematic uncertainties of Union compilation will affect 
strongly the error estimation of dark energy parameters 
[13|. If taking the systematic uncertainties into account, 
we find the limit on the constant EoS of dark energy is 
w = — 0.965±0.080 (1 a) and the error bar is significantly 
enlarged. 

For the time evolving EoS, Wdda) — wq + wi{l — a), 
in Figl5]we illustrate the constraints on the dark energy 
parameters wq and wi. For the flat universe, from the 
latest observational data we find that the best fit model is 



FIG. 2: Constraints on the dark energy EoS parameters wo 
and Wl from the current observations, CMB+BAO+SN. The 
red solid lines and the blue dash-dot lines are obtained for the 
flat and non-flat universe, respectively. And the black dashed 
lines are obtained when (incorrectly) neglecting dark energy 
perturbations. The magenta solid lines stand for uiq = — 1 and 
Wo + Wl = — 1. In this numerical calculation the systematic 
uncertainties of Union compilation is not considered. 



the Quintom dark energy model, wq = —1.08 and wi = 
0.368, whose w{z) can cross the cosmological constant 
boundary during the evolution. But then the variance of 
Wo and wi are still large, namely, the 95% constraints on 
Wo and wi are —1.22 < wq < —0.721 and —1.33 < wi < 
0.947. This result implies that the dynamical dark energy 
models are not excluded and the current data cannot 
distinguish different dark energy models decisively. The 
ACDM model, however is still a good fit right now. 

Our results are consistent with the WMAP5 group 0, 
while the upper limit on wq and lower limit on wi are 
slightly weaker than theirs. This difference is mainly 
from the supernovae datasets used. The supernovae 
dataset we use in this paper is the new Union compilation 
with homogeneous analysis of the present world's super- 



novae data. But in Ref. B th ey use three supcrnovae 
datasets, SNLS [Ml, HST [H and ESSENCE ^]. For 
each of them, they marginahze oyer the absolute magni- 
tude separately, and simply add these three pieces to get 
the total x^- 

Because the parametrization of EoS of the dark en- 
ergy used in this paper is assumed to extend to an arbi- 
trary high redshift, it is important to check if the energy 
density of dark energy component is negligible compared 
with the radiation density at the epoch of the Big Bang 
Nucleosynthesis (BBN), z ^ 10^. As shown by the red 
solid lines in Fig [21 the dynamical dark energy models 
allowed by the current data within the 95% confidence 
level safely satisfy the limits of iuq + wi < and BBN 
|33l . |3J| to avoid the dark energy domination in the early 
universe. 

Furthermore, the Null Energy Condition (NEC) should 
also be satisfied for the EoS of universe w,, (Sal : 



Wu(a) = — =^ — > -1 



(5) 



where Wi and pi are the EoS and energy density for com- 
ponent i in the universe. Violation of NEC will lead to 
the breakdown of causality in general relativit y a. nd the 
violation of the second law of thermodynamics [36| . This 
requirement from Wu{a) > — 1 will constrain the EoS pa- 
rameters of dark energy models [37] . 

Firstly, we consider the WCDM dark energy model. 
From Table I, the current observational data give the 
present energy density of dark energy Jldo — 0.74. If we 
just assume that the EoS of dark energy is constant from 
the early time of universe to present, using the above 
Eq.®, we obtain that the NEC hmit requires the EoS 
of dark energy Wdc ^ —1.35 straightforwardly [3M|. For- 
tunately, the current constraints on WAa with and with- 
out the systematic uncertainties of "Union" compilation 
satisfy this limit safely. However, if we extend the pe- 
riod of validity of constant EoS of dark energy into the 
future where the dark energy component will dominate 
the universe entirely (namely f^dc ~ 1), the NEC requires 
Wdc ^ ^1- Consequently, the phantom dark energy mod- 
els (wac < — 1) will be excluded and the fate of universe 
with Big Rip ^ could not happen. 

For the RunW dark energy model, we illustrate the 
constraints on the dark energy parameters wq and wi 
from the current observations (black solid lines) and the 
NEC (blue dashed lines) in Fig[3l respectively. In the up- 
per panel, the parameterized EoS Wdc — wq + wi{l — a) 
is assumed to be valid until now, not for the future, i.e. 
< a < 1. The reasons for doing this are i) exper- 
imentally there are no constraints on the dark energy 
models for the future universe; ii) theoretically it is al- 
ways possible to choose a different type of parametriza- 
tion of dark energy EoS for the future universe, such 
as Wdc = woCxp(l — a) for a > 1, which matches to 
WriR = Wq + wi{l — a) at a = 1 and satisfies the NEC 
[37|. In fact, if the current EoS of dark energy wq < —1, 
this scenario allows the transition from Wdc < — 1 to 
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FIG. 3: Constraints on the dark energy EoS parameters wo 
and wi from the current observations (black solid lines) and 
the Null Energy Condition (blue dashed lines). In the upper 
panel, the parametrization of EoS in Eq.Q is assumed to be 
a description of dark energy for the past until now. In the 
lower panel, we assume that the parametrized EoS are also 
valid for the far future. 



Wdc > —1, consequently avoids the violation of NEC. 
One can see from the FiglJl in this case, the NEC limit 
does improve the constraints on dark energy parameters 
from the current observations. For example, the current 
observational data permit wq < —1.35 at 2 cr confidence 
level, which however violates the NEC limit. 

In the lower panel of Figl3l the parametrization w^q = 
Wq + wi{l — a) is assumed to be a valid description of 
dark energy model at any time from the far past to the 
far future. One can see from this figure, the NEC puts a 
stronger constraint on the EoS parameters of dark energy 
models than the observational data. For this case, the 
regions of quintessence and phantom dark energy models 
will be shrunk significantly in the {wo,wi) space, and 
the NEC excludes the dark energy models corresponding 
to the regions labelled by the phantom and Quintom A 
with Wdc < — 1 for a — > oo. Consequently, the models 
satisfying the NEC include the quintessence and some of 
the Quintom B dark energy models, which can be seen 
in the lower panel of Fig[31 



TABLE II. Constraints on cosmological parameters rig, aa, 

r, Qk and ^ rUv from the current observations. We have 

shown the mean 1,2 a errors. For the weakly constrained 

parameters we quote the 95% upper hmits instead. 



Parameters 


ACDM 


RunW 


100 X Qk 


A nsi +0-545+1.025 


n nQS+0'8Sl+1.655 
'-'•"yo-o. 881-1. 605 


J2m. 


< 0.533 (95%) 


< 0.974 (95%) 


J2 iTT-u (w/ Lya) 


< 0.161 (95%) 


< 0.252 (95%) 


Us 


n Qfil +0.012+0.024 
U.yDi_g gj2-0.023 


964+0013+0. 027 
^•^'-'^-0.013-0.025 


Cls 


n niQ+0-017+0.032 


n noQ+0. 019+0. 039 


r 


< 0.200 (95%) 


< 0.268 (95%) 




Finally we discuss the degeneracy between the dark en- 
ergy and curvature. As we know, the EoS of dark en ergy 
is degenerated with the curvature of universe fi^ [40|, |4l| . 
If we do not include the prior that the universe is flat, 
from Table I and the blue dash-dot lines in Fig[21 we can 
see that the constraints on wq and wi are weakened sig- 
nificantly and the two-dimensional distribution extends 
more towards the Quintom B region. But the main con- 
clusions are unchanged, wq = —1.11 and wi = 0.475, 
namely, the Quintom dark energy model is still mildly 
favored by the current observational data. Moreover, in 
order to show the importance of dark energy perturba- 
tions in the global analysis, we do the calculation by in- 
correctly neglecting the dark energy perturbations. Il- 
lustrated as the black dashed lines in FigUl one can see 
that the constraints on the dark energy parameters be- 
come tighter immediately, similar to our previous results 
[201 . This study shows that how biased the result will 
be, once the dark energyperturbations are incorrectly 
neglected in the analysis [3, \lM, ^^ ■ 



B. Curvature of Universe 

The measurements on the position of first acoustic 
peak of CMB temperature power spectrum have been 
used to determine the curvature of universe fife. However, 
due to the well-known degeneracy between f2m and flk, 
we have to add other cosmological data, such as the large 
scale structure and supernovae data, to break this de- 
generacy and improve the constraint. Within the ACDM 
model, from Table II and FigUl we can see that our uni- 
verse is very close to flatness, namely, the 95% limit is 
—0.012 < rifc < 0.009, which is consistent with the pre- 
diction of inflation paradigm. 

As we mentioned before, dark energy parameters and 
fife are correlated via the cosmological distance informa- 
tion. In the framework of dynamical dark energy mod- 
els, the constraint on flk should be relaxed. Based on 
the calculations, we can see that the combination of ob- 
servational data implies —0.015 < flk < 0.018 at 2 a 
confidence level. 



FIG. 4: ID current constraints on the inflationary parame- 
ters Us, Qs and r, as well as the curvature flk and the total 
neutrino mass ^ m^, in different dark energy models: ACDM 
model (blue dash-dot lines) and RunW model (red solid lines). 
For ^ m^, we also show the limits combined with the SDSS 
Lyman-Qf forest power spectrum. Blue dotted lines and red 
dashed lines denote the ACDM and RunW models, respec- 
tively. 



C. Neutrino Mass 

Detecting the neutrino mass is one of the challenges 
of modern physics. Currently the neutrino oscillation 
experiments, such as atmospheric neutrinos experiments 
[42] and solar neutrinos experiments [43| , have confirmed 
that the neutrinos are massive, but give no hint on their 
absolute mass scale. Fortunately, cosmological observa- 
tional data can provide the crucial complementary in- 
formation on absolute neutrino masses, because massive 
neutrinos leave imprints on the cosmological observa- 
tions, such as the Hubble diagram, CMB temperature 
power spectrum and LSS matter power spectrum [44| . 

Within the ACDM model, from Table II one can 
read 95% upper limit of the total neutrino mass de- 
rived from the current observations, CMB-I-LSS-I-SN, 
^m,y < 0.533 eV (95% C.L.), which is consistent with 
the recent results from WMAP5 group [SJ, |9|. However, 
there are degeneracies between the neutrino mass and 
other cosmological parameters, such as the EoS parame- 
ters of dark energy [45| and the running of spectral index 
[46| . Due to the degeneracy amon g da rk energy parame- 
ters and the neutrino mass 0, H, [43] , in the framework 
of dynamical dark energy models, the limit on the neu- 
trino mass can be relaxed to 'Y^rrii, < 0.974 eV (95% 
C.L.) significantly, as shown in Figd) 

It is well known that when neutrinos become non- 
relativistic at late time of the universe, they damp the 
perturbations within their free streaming scale. Thus the 
massive neutrinos will suppress the matter power spec- 
trum at small scale by roughly AP/P ~ —8n,,y/V.m [48j . 
Therefore, Lyman-a forest data at small scale can sig- 
nificantly improve the constraint on the neutrino mass. 



But when we use the Lya data, we should keep in mind 
on their unclear systematics right now [^, [2^. When 
including the SDSS Lya forest power spectrum [27l |. 
we can obtain a much more stringent 2 a upper limit 
J2m^ < 0.161 eV in the ACDM model. 

Moreover, the Heidelberg-Moscow (HM) experiment, 
which is controversial for the time being, uses the half- 
life of Oz^2/3 decay to constrain the effective Majorana 
mass and this translates to the constraint on the sum of 
neutrino masses under some assumptions [49[ , ^ m^ ~ 
1.8 ±0.6 eV (95% C.L.). We can find an obvious tension 
on the neutrino mass limits from between the cosmologi- 
cal observations and the HM experiment, which however 
can be resolved if the neutrino masses vary during the 
evolution of the universe [50| . In order to be consistent 
with the observational data, the neutrino mass must be 
very small in the past, but has grown recently in order 
to agree with the HM experiment data. 

Again, in the RunW model, the limit on the neutrino 
mass is relaxed, ^ rrii, < 0.252 eV (95% C.L.), which im- 
plies the existence of degeneracy between the dark energy 
parameters and neutrino mass. 




FIG. 5: 68% and 95% constraints on the panel {ns,r) based on 
the different Dark Energy models: ACDM model (blue dash- 
dot lines) and RunW model (red solid line). The two magenta 
solid lines delimit the three classes of inflationary models, 
namely, small-field, large-field and hybrid models. The star 
points are predicted by HZ spectrum, m^(/)^ model and A0* 
model, respectively. These predictions assume that the num- 
ber of e-foldings. A'', is 50 — 60 for rn?4>'^ model and 64 for A</f)* 
model. 



D. Inflationary Parameters 



Inflation, the most attractive paradigm in the very 
early universe, has successfully resolved many problems 
existing in hot Big Bang cosmology, such as flatness, hori- 
zon, monopole problem and so forth 51*. Its quantum 
fluctuations turn out to be the primordial density fluc- 
tuations which seed the observed large scale structures 
and the anisotropics of CMB. Inflation theory has suc- 
cessfully passed several non-trivial tests. Currently, the 
cosmological observational data are in good agreement 
with a gaussian, adiabatic and scale-invariant primordial 
spectrum, which is consistent with single field slow-roll 
inflation predictions. 

Within the ACDM model, from Figdl we obtain the 
limit on the spectral index of Ug — 0.961 ± 0.012 (1 cr), 
which excludes the scale- invariant spectrum, ris = 1, and 
the spectra with blue tilt, Us > 1, at more than 3 a con- 
fidence level. When considering the gravitational waves, 
the latest observational data yield the 95% upper limit 
of tensor-to-scalar ratio r < 0.200. In FigH] we show the 
two dimensional constraints in (ris ,r) panel which can be 
compared with the prediction of the inflation models. We 
find that the Harrison-Zel'dovich-Peebles scale-invariant 
(HZ) spectrum {ug = 1, r = 0) is still disfavored more 
than 2 a confidence level. And many hybrid inflation 
models and the inflation models with "blue" tilt (n^ > 1) 
are also excluded by the current observations. Further- 
more, the single slow-rolling scalar field with potential 



¥{<!>) 



which predicts {ns,r) = (1 - 2/N,8/N), 



is still well within 2 a region, while another single slow- 
rolling scalar field with potential F(0) ~ At/)^, which pre- 
dicts {ns,r) = (1 — 3/A^, 16/A^), has been excluded more 
than 2 cr fifll,!!!. 



However, the tensor fluctuations and the dark energy 
component, through the ISW effect, are correlated, which 
mostly affect the large scale (low multipoles) temperature 
power spectrum of CMB [l|, [53| . In the framework of dy- 
namical dark energy model, we find that the upper limit 
of r can be relaxed to r < 0.268 (95% C.L.). Further- 
more, the 95% confidence level contour in (ns,r) panel 
will be enlarged consequently and the distribution ex- 
tends towards the hybrid inflation region. Therefore, we 
can see that the HZ spectrum is consistent with the latest 
observational data and many hybrid inflationary models, 
the inflationary models with "blue" tilt {ug > 1), which 
are excluded in the ACDM model, have revived in the 
framework of dynamical dark energy model as illustrated 
in FigO 

Finally, we explore the constraint on the running of 
the spectral index from the latest observational data. 
When combining the WMAPl or WMAP3 data with 
other astronomical data, the previous analysis have found 
a significant evidence for large running |54l . l5q . Physi- 
cally this large running would be a great challenge to the 
single- field inflation models |56l . |57| . The combination of 
the current observational data yield the limit on the run- 
ning of the spectral index of ag = —0.019 ± 0.017 (1 a) 
for the ACDM model and a^ = -0.023 ± 0.019 (1 a) for 
the RunW model, respectively. The error is dramatically 
reduced compared with the previous results [l|, |53| , ben- 
eficial from the more accurate observational data. Given 
the current data, we find no significant evidence for the 
large running of the spectral index. 



IV. SUMMARY 

Recently many experimental groups have published 
their new observational data, such as temperature and 
polarization power spectra of WMAP5 [1, [3, Uol, tem- 
perature power spectrum of ACBAR [ll| and supernovae 
dataset of Union compilation [l3| . In this paper we report 
the updated constraints on the cosmological parameters 
from these latest observational data, such as the EoS of 
dark energy, curvature of universe, neutrino mass and 
inflation parameters. 

For dark energy, we explore the constraints on two 
kinds of dynamical models. Assuming a flat universe, 
within the WCDM model, we find that the latest obser- 
vational data yield the limit on the constant EoS of dark 
energy, w = -0.977 ± 0.056 (1 a). For RunW model 
with a flat universe, we find that the best fit model is 
wq = —1.08 and wi — 0.368 described by the Quintom 
theory with EoS across the cosmological constant bound- 
ary. And because the precision of current observations 
are not good enough to determine the dark energy EoS 
conclusively, the dynamical dark energy models are not 
excluded and the ACDM model remains a good fit. 

The Quintom scenario, with the particular feature that 
its EoS can cross the cosmological constant boundary 
smoothly, has been applied to many aspects of cosmology 
theoretically. Firstly, the Quintom dark energy models 
can also give rise to interesting prediction on the fates of 
the universe, different from the Quintessence or Phantom 
models, such as the cyclic universe [5^^], an expanding 
universe with oscillating EoS. Secondly, applying a Quin- 
tom matter for the early universe can provide a scenario 
of bouncing cosmology, which can avoid the notorious 
issue of initial singularity [60| . 

Our results also show that the universe is very close to 
flatness and the upper limit on the total neutrino mass 
is X^'^i' < 0.533 eV (95% C.L.), from the combina- 
tion of CMB, BAO and SN data. Given the efficiency 
of Lya forest data on constraining the total neutrino 
mass, we also perform a calculation with the inclusion of 
the SDSS Lyman-a forest power spectrum and find that 



J2m^ < 0.161 eV (95% C.L.). This result might lead 
to the exclusion of degenerate pattern of neutrino mass, 
when combining the results of neutrino oscillation experi- 
ments. Due to the degeneracy between the neutrino mass 
and EoS of dark energy, however in the presence of dy- 
namics of dark energy, the constraints on ^ rrii, can be 
relaxed by a factor of 2. 

Finally, for the inflationary models, within the ACDM 
framework, we find that the latest observational data pre- 
fer the inflation models with "red" tilt, namely, Ug = 
0.961 ± 0.012 (1 a) and small tensor fluctuations, r < 
0.200 (95% C.L.). Because of the degeneracy between r 
and EoS of dark energy, the upper limit on r is relaxed to 
r < 0.268 (95% C.L.) and the parameter space in {ns,r) 
panel are enlarged in the framework of dynamical dark 
energy models. Therefore, the inflationary model with 
HZ primordial spectrum (rig = 1, r = 0), some hybrid 
models and some models with a "blue" tilt {ug > 1), 
which are excluded more than 2 a confidence level in the 
ACDM model, will be consistent with the current ob- 
servations now. Furthermore, in our analysis we do not 
find any significant evidence for the running of spectrum 
index. 
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